Tunable bandgaps and excitons in doped semiconducting carbon nanotubes made 

possible by acoustic plasmons 
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Doping of semiconductors is essential in modern electronic and photonic devices. While doping 
is well understood in bulk semiconductors, the advent of carbon nanotubes and nanowires for na- 
noelectronic and nanophotonic applications raises some key questions about the role and impact of 
doping at low dimensionality. Here we show that for semiconducting carbon nanotubes, bandgaps 
and exciton binding energies can be dramatically reduced upon experimentally relevant doping, and 
can be tuned gradually over a broad range of energies in contrast to higher dimensional systems. 
The later feature is made possible by a novel mechanism involving strong dynamical screening effects 
mediated by acoustic plasmons. 



Nanomaterials have been lauded for their promise in 
electronic and photonic applications [l|. Quite often, 
the imagined nanodevices rely on analogies with those 
based on bulk semiconductors. However, the true po- 
tential of nanomaterials lies in the exploitation of their 
unique properties to realize entirely new device concepts. 
In particular, approaches for externally controlling their 
electronic and optical properties would enable new strate- 
gies for device design. 

Here we propose that such control is possible in car- 
bon nanotubes (CNTs) through electrostatic doping. We 
find that quasiparticle (QP) band gaps and exciton bind- 
ing energies can be reduced dramatically by hundreds of 
meVs upon doping, and yet prominent optical absorption 
features shift by relatively small amounts. Furthermore, 
we show that doping has a unique influence on CNT ex- 
citon properties: in contrast to bulk excitons, bound ex- 
citons in semiconducting CNTs are not quenched by dop- 
ing and their binding energy can be tuned gradually even 
at very high doping. These features arise due to the pres- 
ence of acoustic plasmons and their impact on dynamical 
screening. 

We utilize a many-body ab initio approach [24j] to 
calculate the electronic and optical properties of electro- 
statically doped semiconducting CNTs. We focus on the 
semiconducting (10,0) CNT, with diameter D = 0.78 nm, 
and perform ab initio calculations [5| at zero doping and 
for a free carrier concentration p = 0.6 holes/nm. We 
use the GW approach [2] to obtain QP properties near 
the r point and solve the Bethe-Salpeter (BS) equation 
for the excitonic effects [3||- Applying this approach to 
doped CNTs necessitates careful consideration because 
of the presence of acoustic plasmons, a unique feature of 
low-dimensionality materials [g, |7| . 

Indeed, in quasi-lD systems such as CNTs, electron 
gas and tight-binding models predict "acoustic" plas- 
mons, whose energies approach zero in the long wave- 
length limit Ldapiq — ^ 0) oc q^/p\og{\q\D/2). Our ab 
initio calculations also reveal these plasmons in doped 
CNTs [8(: Fig. [1^ shows the inverse dielectric func- 
tion Eof^iq = 0.35nm-\£') of the (10,0) CNT at p = 



0.6 holes/nm. The peak in Ime~^ signals a low-energy 
plasmon, which gives rise to a transition in Ree~^ be- 
tween a very small value at zero energy and a value 
close to 1 above the plasmon energy, i.e. a transition 
between metallic-like and semiconducting-like screening 
9] . These acoustic plasmons span a broad range of ener- 
gies, as seen in Fig. [TJd, and dynamical screening effects 
due to them are very important for both QPs and exci- 
tons and cannot be neglected or simply integrated out. 

To include these dynamical screening effects during the 
GW calculations we modify methods previously applied 
in the context of doped bulk semiconductors [lOf . In the 
undoped case, optical plasmons in CNTs have energies 
above '^ 5 eV ll|, and their contributions is taken into 
account through the Generalized Plasmon Pole (GPP) 
approximation [2]. In the doped case, GPP alone does 
not describe satisfactorily dynamical effects from acous- 
tic plasmons and the screening is evaluated instead from 
e~^ = e~^^ + Se~^, where Smt is the dielectric function 
of the intrinsic semiconductor for which we make use of 
the GPP approximation and Se^^ is obtained within the 
Random Phase Approximation. 

Dynamical effects are included in the BS calculations 
by an effective dielectric screening i that depends self- 
consistently on the binding energy Eb of the exciton [l^l . 
Taking advantage of the fact that for the excitons consid- 
ered here most of the corresponding electron- hole transi- 
tions have energies close to the onset of the electron-hole 
continuum, and neglecting finite lifetime effects, one can 
write to first order in dynamical effects [13| : 
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When low-energy plasmons are absent, the commonly 
used static approximation is obtained by retaining the 
first term on the r.h.s. of Eq. ([T]). The second term 
captures dynamical effects due to acoustic plasmons (we 
neglect dynamical effects due to optical plasmons during 
the BS calculations in both undoped and doped cases). 
Smaller, higher order corrections in dynamical effects 
(not shown in Eq. ([T])) due to acoustic plasmons are 
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FIG. 1: Dielectric properties of the (10,0) doped CNT. a) 
Real and imaginary parts of the inverse dielectric function 
as a function of energy calculated ab initio for a fixed wave- 
vector q = 0.35 nm^^. The subscripts on e^^ imply their 
evaluation at G = G' = 0. b) Plasmon dispersion relation 
calculated ab initio and within the model screening used in 
the effective mass model. Doping level p — 0.6 holes/nm. 



included as well in our calculations. 

Fig. [5] shows a sketch of the calculated QP bands and 
exciton level associated with the lowest prominent opti- 
cal transition (£'ii), before and after doping, where the 
energy scale has been preserved between the two cases. 
The quantities of interest are the bandgap E^^ and the 
exciton binding energy E]^. It is clear from this fig- 
ure that the bandgaps and exciton binding energies are 
significantly reduced by doping. In fact, at this doping 
(p = 0.6 holes/nm), Eg^ is reduced by 800 meV, while 
E^ is reduced by 590 meV. Both of these changes are 
large by any measure; in particular we estimate that band 
gap renormalization (BGR) is about an order of magni- 
tude larger than in typical bulk semiconductors at the 
same doping [14] . The QP bands and exciton level associ- 
ated with the second lowest prominent optical transition 
{E22, not shown in Fig. ^ also suffer from a significant 
bandgap reduction of 130 meV '15^, and a decrease of 
the exciton binding energy by 240 meV. Moreover, dop- 
ing leads to an important 40% reduction in the effective 
mass at the valence band maximum for En. 

The excitonic effects are presented in more detail in 
Fig. 131 where the optical absorption spectra for light po- 
larized along the tube axis is calculated in two ways. The 
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FIG. 2: Illustration of bandgap and exciton renormalization 
in the (10,0) CNT. The blue lines indicate the En valence 
and conduction bands, and the bound exciton is indicated 
by a red line. Occupied electronic states are indicated by 
the shaded areas, and the highest occupied energy levels have 
been aligned between the doped and undoped cases. The dif- 
ferent band curvatures represent the change in effective mass 
upon doping. 



blue lines are obtained neglecting the electron-hole inter- 
action, with the corresponding onset energies indicating 
the electron-hole continua. The red lines show the cor- 
rect absorption spectra obtained with the electron-hole 
interaction included. As each bound exciton requires a 
separate self-consistent calculation in the doped case, we 
have focused on the lowest energy bright exciton associ- 
ated with each of the En and E22 continua. 

In the undoped case, the En and E22 excitons show 
very large binding energies: E]^{p = 0) = 0.79 eV 
and E'^{p = 0) = 1.00 eV. As discussed above, upon 
doping, a dramatic change in excitonic properties oc- 
curs. As seen in Fig. [3}d, E]^ suffers a decrease of 
- 0.6 eV to E}j^{p = 0.6 holes/nm) = 0.20 eV, while 
the corresponding E^ renormalization is 0.24 eV (see 
Fig. [3ji). While both the En and E22 excitons are 
affected by the change in dielectric screening, the En 
exciton renormalizes more because it is affected by the 
bleaching of transitions. We also note from Fig. [3}d a 
six- fold reduction in the oscillator strength of the En 
exciton, in good accord with recent photo-luminescence 
measurements [l6| which assigned a factor of five in the 
drop of the En exciton radiative decay rate of a 1.4 
nm diameter CNT [for an estimated maximum doping 



Pmax ~ 0-16 holes (electrons) /nm]. 

Our ab initio results suggest large changes in bandgap 
and exciton properties upon doping 17[. This could lead, 
for example, to engineering of CNT optoelectronic de- 
vices by electrostatic control; but taking advantage of 
these new features requires robust control of doping- 
induced properties. Because our many-body ab initio 
calculations are extremely demanding, exploring the tai- 
loring over a broad range of doping is not possible. There- 
fore, we developed a compact model for excitons in doped 
CNTs based on an effective mass (EM) approach [18|. 

Our EM model seeks to describe the binding energy 
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FIG. 3: Impact of doping on the electronic and optical prop- 
erties of the (10,0) CNT. a,b) Optical absorption spectra near 
the Ell transition with (red) and without (blue) the electron- 
hole interaction, a) Before doping, b) After doping with 0.6 
holes/nm. c,d) Optical absorption spectra near the E22 tran- 
sition c) Before doping, d) After doping with 0.6 holes/nm. 
Arrows indicate exciton binding energies. Optical spectra are 
broadened by a 20 meV Lorentzian. 



and envelope function of excitons in CNTs taking ad- 
vantage of their large spatial extent along the tube axis 
relative to inter-atomic distances. The effective interac- 
tion between the electron and hole composing an exciton 
depends on e via [see Eq. ([T|)] the dielectric function 
e = 1 — vP, where P is the irreducible polarizability. 
Our EM approximation assumes that P{z, z', p, p', tp, ip') 
is localized on a cylindrical tubule with radius R, and 
that local field effects along the tube axis can be ne- 
glected. Noting that only the L=0 component of the 
angular momentum L is relevant for screening effects in 
excitons composed of electron-hole transitions between 
bands of same L, one replaces P{z,z' , p, p' ^L = 0) with 
P{z — z',R,L — 0). Relevant dynamic effects are in- 
cluded via P{uj) « Pf'''"'{uj) + P™*(a; = 0), where 
pfree (jgscribes intraband transitions within parabolic 
band approximation [6|, and P™* represents the inter- 
band transitions for the intrinsic semiconductor. (We 
keep transitions from the top of valence band to con- 
duction states as wc found that they make very small 
contributions.) p*"* is extracted from ab initio calcu- 
lations for the intrinsic polarizability (p™*) by impos- 
ing that P*"* and P*"* yield the same average along 
the radial direction over one unit cell. In Fourier space 

this reads: Ptloi^l) = To ^S^g' =0^"^^ ' ""^^'''^ ^"'^ '^ ^^^ 
cross-sectional area of the unit cell considered in the ab 
initio case. The ability of our model -free of any ad- 
justable parameters- in describing dynamical screening 
effects from acoustic plasmons is demonstrated in Fig. 
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tended to carrier densities as small as p « 1 hole/800 nm 
(in terms of number of holes per atom, this corresponds 
to 10^^ — 10^® holes/cm in typical bulk semiconductors, 
i.e., approximately the degenerate limit). Comparison 
with ab initio results illustrates the high accuracy of our 
EM model. Fig. \^ shows that with good approxima- 
tion, SEf = Ef{p) - Sf (0) ex ^, which we found 
to be a signature of acoustic plasmons. A similar trend 
is found for 5E^ , where deviations from the y/p behav- 
ior are more pronounced due to bleaching of transitions. 
More importantly, the mild dependence on doping should 
be contrasted to that in higher dimensional semiconduc- 
tors where, at equivalent doping levels, excitons are either 
quenched [19| or the dependence on doping is exponential 
[20|, giving poor control over optical properties. 

The large change in exciton binding energy combined 
with the mild dependence on doping implies that exci- 
tonic properties in CNTs can be efficiently controlled 
through doping. The origin of this feature lies in the 
presence of acoustic plasmons. To emphasize this point. 
Fig. |3b also shows EM results for exciton binding ener- 
gies within the static approximation, i.e. without proper 
inclusion of dynamical effects due to acoustic plasmons. 
In this case, binding energies drop exponentially with 
doping, much like is observed in two-dimensional ma- 
terials [20[; moreover, the En exciton gets quenched 
beyond p « lhole/15 nm. The importance of acous- 
tic plasmons is readily seen within our model and from 
Eq. (dJ, which can be shown to yield: 
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For long 



wavelengths where uiapil) << Eb, one has e{q;EB) ~ 
£mt('Z,0), as opposed to the static approximation result 
e{q; Eb) ~ s{q, 0). The difference between these two val- 
ues can be orders of magnitude depending on q. 

We obtain the QP fundamental bandgap versus dop- 
ing using our exciton EM model and our many-body 
ab initio results for optical properties. Indeed, we can 



write E^^ = 






f7" - AE},^ where fi" is the ex- 

^cont ^g 

Fig. [2])- We calculate the doping dependence of AE^p 



citon excitation energy and AE^ 



(see 



from AE]} = 



kl/2p* 



with the reduced exciton mass 



p* obtained at various doping levels by interpolating the 
values from many-body ab initio calculations at p = 
and p — 0.6 holes/nm. Furthermore, we find from our 
ab initio calculations at p = 0.6 holes/nm that ^^^ only 
increases by ~ 0.1 eV upon doping (similarly for fi^^), a 
result of cancelation between self-energy corrections and 



Fig. S^ shows the EM results for E]^ and E"^ ex- 



excitonic effects [2l| . The smallness of these shifts is in 
excellent agreement with measurements of the £'33 ab- 
sorption peak of a 1.4 nm diameter CNT jl6[, showing 
a red-shift of 20 meV at p^ax- Therefore, with little ex- 
pected error, we assume a linear dependence of f2^^ on 
doping, and plot in Fig. Hh the fundamental bandgap 
versus doping; the trend indicates that BGR can also be 
tuned gradually over a broad energy range. 
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FIG. 4: Exciton binding energies and fundamental bandgap 
as a function of doping in the (10,0) CNT. a) Binding en- 
ergies for the Ell and E22 lowest bright excitons calculated 
from the EM model and compared to ab initio results, b) 
Change in exciton binding energies, SE% = Eg{p) — Eg{0), 
as a function of \J pj ptot where ptot is the total electron charge 
density, within the EM model with and without dynamical ef- 
fects. The dotted lines refer to the EM model. The triangles 
represent ab initio results: within the static approximation, 
we find no bound exciton associated with En and E"^ ~ 30 
meV at p = 0.6 holes/nm. c) Dependence of the fundamental 
bandgap on doping, based on the EM model and compared 
to ah initio results. 



The giant BGR discussed here is in agreement with 
recent experimental results [22| that combined photocur- 
rent spectroscopy with transport measurements: for a 
1.5 nm diameter CNT with QP band gap of 0.91 eV at 
zero doping, a BGR of 0.54 eV was deduced at a doping 
density of 0.7 electrons/nm. 

In conclusion, we have shown that doping has a pro- 
found and unique impact on the electronic and optical 
properties of semiconducting CNTs, and that dynamical 
effects from acoustic plasmons are essential to capture 
these effects. Our study indicates that bandgaps and 
exciton binding energies in CNTs can be tuned signifi- 
cantly and gradually by electrostatic doping, establish- 
ing a new framework for the understanding and design 



of CNT-based devices. We expect that similar control 
will be possible in a broad range of nanomaterials. 
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